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Abstract 

Light transmission or diffraction from different quantum phases of cold atoms in an optical 
lattice have recently come up as a useful tool to probe such cold atomic systems. The periodic 
nature of the optical lattice potential closely resembles the structure of a diffraction grating in 
^ \ real space, but loaded with a strongly correlated quantum many body state which interacts with 



the incident electromagnetic wave, a feature that controls the nature of the light transmission 
or dispersion through such quantum medium. The paper contains a detailed analysis of such a 



quantum diffraction grating. We have particularly studied how the dispersion shift induced by the 



cavity changes as a function of the relative angle between the cavity mode and the optical lattice. 
PACS numbers: 03.75.Lm,42.50.-p,37.10.Jk 
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I. INTRODUCTION 

Ultra cold atomic condensates loaded in an optical lattice provide a unique op- 

portunity to study the properties of an ideal quantum many body system. After the first 
successful experiment in this field js] where a quantum phase transition from a Mott Insula- 
tor(MI) to Superfiuid(SF) phase was observed, extensive theoretical as well as experimental 
study in this direction took place. The field continues to be a frontier research area of atomic 
and molecular physics, quantum condensed matter systems as well as quantum optics si- 
multaneously, and holds promise for application in fields like quantum metrology, quantum 
computation and quantum information processing 0, 5|. 

The relevance of the field of ultra cold atomic condensates to quantum optics was sug- 
gested much earlier when it was pointed out that the refractive index of a degenerate Bose 
gas gives a strong indication of quantum statistical effects 6| and the interaction between 
quantized modes of light and such ultra-cold atomic quantum many body system is going 
to lead to a new type of quantum optics j7|. Subsequently, it was pointed out that the 
optical transmission spectrum of a Fabry-Perot cavity loaded with ultra cold atomic BEC 
in an optical lattice can clearly distinguish between a SF and MI state Q and may be used 
as an alternative way of detecting such phase transition without directly perturbing the 
cold atomic ensemble through absorption spectroscopy. A successful culmination of some 
of these theoretical predictions happened with the recent experimental success of realiza- 
tion of a strongly coupled atom-photon system where an ultra cold atomic system is placed 
inside an ultra-high finesse optical cavity jo], [lo|, such that a photon in a given quantum 
state can interact with a large collection of atoms in same quantum mechanical state and 
thereby enhancing the atom-photon coupling str eng th. As an aftermath, a host of inter- 
esting phenomena such as cavity optomechanics observation of optical bistability and 



Kerr nonlinearity 



12| has been experimentally achieved with such systems. It may be also 



mentioned in this context that Bragg diffraction pattern from a quasi-two dimensional Mott 
Insulator, but without any cavity, was also recently been observed experimentally 13|. 

The theoretical progress in understanding such atom-photon systems involving ultra cok 
atomic condensates is also impressive. A series of work by the Innsbruck group 
clearly pointed out how the optical properties of the cavity reveals the quantum statistics 
of these many body systems. In another set of work, cavity induced bistability in the MI 
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to SF transition either due to strong cavity-atom coupling [2l| or due to the change in the 



boundary condition of the cavity 
optomechanics 



23 



22| has been studied and its relation to cavity quantum 



24j has also been explored. The self organization of atoms in a mul- 



timode cavity due to atom-photon interaction leading to the formation of exotic quantum 
phases and phase transition 
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27j is another major development in this direction. The 



recent observation of Dicke quantum phase transition through which a transition to a super- 



28| through such self organization is an important experimental 



sohd phase was achieved 
landmark in this direction. 

The physics of ultracold atoms loaded inside an ultrahigh-finesse Fabry- Perrot cavity can 
be analyzed from two different, but highly correlated perspectives. For example, ultra cold 
atomic ensemble loaded in such optical lattices with short range interaction can exist in two 
different types of quantum phases, MI and SF. The former is a definite state in the Fock space 
with well defined number of particles at each lattice site and lacks phase coherence between 
the atomic wavefunctions at different sites. The latter is a superposition of various Fock space 
states and has phase coherence. A phase transition between these phases takes place as the 
lattice depth varies. The statistical distribution of atomic number that characterizes these 
many body states consequently influences the transmitted or diffracted electromagnetic wave 
through atom-light interaction and thereby changes the dielectric response of such cavity in 
the same way as the change of material leads to the change in refractive index. 

From another perspective, the periodic optical lattice potential forms a grating like struc- 
ture in the real space, but now each slit of the grating contains ultra cold atoms in their 
quantum many-body state, that interacts with the light quanta of the electromagnetic field 
through the dipole interaction. Such a system has been dubbed as a quantum- diffract ion 
grating in the literature It is well known that any quantum mechanical scattering process 
leads to the diffraction effect and thus such effect is ubiquitous in various quantum systems. 
As early as in 1977 in a review article by Frahn 29| an overview of such wide range of 
quantum mechanical diffraction process was presented in a common theoretical framework 
by comparing them with classical optical diffraction. Though some element of such quantum 
diffraction is also present in atom-photon system under consideration, it is unique in the 
sense that here electromagnetic wave is getting diffracted by a quantum phase of matter 
wave loaded inside a cavity. A classical description of such diffraction of electromagnetic 
wave by a single atom or an atomic ensemble placed inside a cavity was also discussed in 
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detail in ref. 



3Q|. 



It has been pointed out 



3l| that in ultra cold fermionic atoms, the diffraction properties 



of such quantum-diffraction grating is strongly dependent on the mode of quantization of the 
cavity. In the limit of very large cavity detuning also, the features of such quantum diffraction 
: or ultra cold bosonic atoms and its departure from the classical behavior has been studied 
1^ . The result from these earlier studies indicate that a detailed analysis of the diffraction 
properties of such quantum diffraction grating has the potential to characterize the many- 
body quantum states of ultra cold atoms in more detail. Since the relevant experimental 
system is already available, such a study is even more encouraging. In the current paper we 
carry an analysis of the diffraction property of such quantum diffraction grating in detail. 

The plan of the paper is as follows. In the next section we begin with a brief review of the 
formalism that is used to calculate the transmission from such a cavity loaded with ultra cold 
atoms in optical lattice. We shall particularly discuss the cases when a single mode and two 
modes are excited in the Fabry-Perot cavity and show how the transmitted intensity through 
the cavity will be calculated in these two cases. The particular emphasis in this work will be 
on how the transmitted intensity changes as one varies the relative angle between the optical 
lattice and pump laser(s) that excite the modes. Using this theoretical framework in the next 
sections we shall present our result which will show how the transmitted intensity varies both 
for the MI and SF phases under different conditions. An analysis of these results and their 
comparison against various classical diffraction pattern provides us a sound understanding of 
this quantum diffraction phenomena. Finally we extend the results to a ring shaped cavity 
and will show how cavity quantization procedure changes the transmitted intensity. 



II. MODEL 



The physical system we describe is depicted in Fig. [Hand consists of identical two- level 
bosonic atoms placed in an optical lattice of M sites inside a Fabry-Perot cavity. K sites 
among these are illuminated by cavity modes, pumped into the system by external lasers. 
We shall consider both the cases, in which single cavity mode, and two cavity modes will be 
excited. 
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FIG. 1: (Color Online) Schematic diagram for cold atoms in an optical lattice loaded in a cavity. 
The optical lattice is created from two counter propagating laser beams and has a site spacing of 
length d. The two standing wave cavity modes, MODE 1 and MODE 2 are at angles 9q and 9i 
with the axis of the optical lattice respectively. The MODE 1 is being pumped by a pump laser 
with amplitude ijq while MODE 2 is not being pumped and is used to collect the scattered photons 
by a detector. In the single mode case, the detector is also on MODE 1 and has not been shown 
in this figure. 



These modes can be composed of either standing waves(SW) or traveUing waves(TW). 
The setup given in Fig. [1] can reahze standing waves; whereas later we shall discuss the 
corresponding setup for travelling wave solutions. Such a system was studied in 
and for a detailed treatment, the reader may refer to [14^]. Here we describe this theoretical 
framework briefly. 

The above mentioned system can be theoretically modeled as a collection of two level 
atoms, that are approximated as linear dipoles, to account for their interaction with the 
quantized electric field of the cavity modes. To describe the system through an effective 
Hamiltonian one then uses well known rotating wave approximation in which all fast os- 
cillating terms in the Hamiltonian are neglected. The excited state of the two level atoms 
is then adiabatically eliminated assuming that the cavity modes are largely off resonant to 
energy difference between the atomic levels. Thus in the resultant system all the atoms are 
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in their ground states. The effective Hamiltonian arrived in this way is given by 



\^ ( K 

a — 



+ E ^ ( E ) + f E %(% - 1) 



(1) 



Here 



= '^huiajai - ih^{r]*{t)ai - rii{t)a\) 



where first term denotes the free field Hamiltonian and the second depicts the interaction of 
classical pump field with cavity mode, ai is the annihilation operators of light modes with 
the frequencies w;, wave vectors k^, and mode functions ui{t). T]i{t) = 7796"*'^''* is the time 
dependent amphtude of the external pump laser of frequency Up that populates the cavity 
mode. 

Here Jj^f. correspond to the matrix element of the atomic Hamiltonian in the site localized 
Wannier basis, w{r — Tj), namely 

= J dvwiv - r,)Hawir - r,) (2) 

where Ha = — + Vci{r) is the Hamiltonian of a free atom of mass rria in an optical lattice 
potential Vc/(r). Therefore, Jq = Jf^, and J'^' =Jf,j^i are respectively the onsite energy and 
the hopping amplitude of proto-type Bose Hubbard model given in ij. At the atomic site 
j, hj is the annihilation operator, and hj = P-hj is the corresponding atom number operator, 
N = J^'jLi % denotes the total atom number and B = J2jLi 

The coefficients Jj™ is similar to J|'^, but now generated from interaction between atoms 
and quantized cavity modes and is given by , 

Jj;^ = J drw{r - r,)u*i{r)um{r)w{r - r^) (3) 

Ala = - denotes the cavity atom detunings where Ua is the frequency corresponding to 
the energy level separation of the two- level atom and g is the atom-light coupling constant. 
Thus the fourth and fifth term in ([1]) respectively contribute to the onsite energy and hopping 
amplitude due to the interaction between atoms and quantized cavity modes. 
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In the last term, U = ^^^^^ J (ir|w(r)|^ , where denotes the s-wave scattering length 
and gives the onsite interaction energy. For a sufficiently deep optical lattice potential Vci{r), 
the overlap between Wannier functions can be neglected. In this limit J'^' = and Jj*^ = 
for j 7^ k. Such Wannier functions can be well approximated as delta functions centered at 
lattice sites and consequently Jj™ = Ui{rj)um{rj). 

The above Hamiltonian in ([T]) describes the zero temperature quantum phase diagram of 
ultra cold bosonic atoms loaded in an optical lattice placed inside a optical cavity. This is 
because their many body quantum mechanical ground state can exist in various quantum 
phases jl, 3] as a function of parameters like U and J. In the subsequent analysis in this 
work, the physical system that diffracts the photons is therefore a novel type of quantum 
diffraction grating not only because the diffracting medium corresponds to a quantum phase 
of ultra cold atoms, but also due to the fact that it is embedded in a optical lattice/grating 
like structure in real space which in turn affects the nature of such quantum phase. As we 
shall point out, one particular way of understanding the nature of such quantum diffraction 
and differentiate from classical diffraction or any other quantum diffraction 29|] is to study it 
as a function of the relative angle between the cavity mode and direction of the optical lattice 
in which the cold atoms are loaded. Quantum diffraction of electromagnetic wave by such 
ultra cold atomic condensate inside an Fabry-Perot cavity, but without loading them in an 
optical lattice (other than the one dynamically generated due to cavity-atom coupling), was 
already experimentally studied in ref ll| in the context of cavity quantum optomechanics. 
Thus the physical system under consideration is very much realizable experimentally. We 
start our discussion by briefly outlining the relevant theoretical framework to understand 



such quantum diffraction following ref 



li. 



III. METHODOLOGY 

From the atom-photon Hamiltonian ([T]), the Heisenberg equation of motion of the photon 
annihilation operator ai is given by 

ai = -iujiai - iSiDiiCii - iSmDimCim - i^ai + r]i{t) (4) 

with Dim = '"«*(rj)wm(rj)nj, where / 7^ m and 61 = g'^/Aia- k is the cavity relaxation 

rate introduced phenomologically. The ffist, fourth and the fifth terms on the right hand side 
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correspond to property of light transmission through an empty cavity. The second and third 
terms give the information about the atom-hght interaction in the cold atomic condensates. 



A. Single Mode 

First we shall consider the case when a single cavity mode is excited. From the stationary 
solution for one mode case, namely ai — we obtain the expression for the corresponding 
photon number operator as, 

dido = (5) 

Here Ap = uop — cuq is the probe-cavity detuning and di = do. In this case dm = and 
Di,m — Dqq. The single mode transmission through the cavity is calculated by taking the 
expectation value of the above expression in given many-body atomic ground state. As 
expected such an expression is similar to the standard Breit Wigner form. However Dqq is 
in terms of the Fock space operators acting on the atomic ensemble, revealing the statistical 
properties of the quantum matter of ultra cold atoms. 

Now in the denominator of the above expression the shift in frequency is determined by 
the eigenvalue of the operator Dqo, which is dependent on both the atomic configuration, i.e., 
the number of illuminated atoms and the mode functions. For plane standing waves the mode 
function, u{rj)sw — cos(k.rj + 0) where is constant phase factor which has been set to 
zero and denotes the position vector of the jth site on the optical lattice. Here we consider 
a one dimensional optical lattice with site spacing d. For a cavity mode of wavelength A 
incident at an angle 6 with the optical lattice, u{rj)sw = cos{^jdcos9). We assume the 
cavity mode wavelength to be 2d and thus, the mode function is u{rj)sw — cos{j7rcos9), 
where j e /. For such standing waves, the factor Doq becomes, 

■Dqo = uluiUj = cos^{j7rcos9)nj (6) 

j j=l:K 

where K are the number of illuminated sites. This shows that shift in the cavity resonant 
frequency is dependent on the relative angle of the cavity mode with the optical lattice. 
To simphfy the analysis, here it has been assumed that while changing this angle, the light 
beam waist is modified in a way that we always illuminate only fixed K sites. However, as 
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explained below, a few sites fall at intensity minima of the cavity mode, thus changing the 
effective number of illuminated sites. 

For example in Fig. [2a] we show the cases when light with wavelength(=2(i) is incident 
at ^ = 0° and 9 = 60°. When the angle ^ = 0°, all the atoms are at the points of 
maximum intensity or the anti-nodes of the cavity mode wavelength. Thus all the atoms 
are illuminated. When 6 changes to 60°, the projected wavelength along the optical lattice 
direction changes and a few atoms which were at the maxima points are now placed at 
the points of minimum(or zero) intensity or nodes. Thus now only the alternate sites are 
illuminated as can be seen in Fig. [2al Therefore, the effective number of illuminated sites 
in the lattice at 6' = 60° reduces to half its value at = 0°. Hence the dispersive shift varies 
with the change in the relative angle of the cavity mode and the optical lattice. 



1. Mott Insulator 



We shall first consider the case when the ground state of the atomic ensemble is a MI, 
ie., a single state in the Fock space. 

= \n, n, n, n) (7) 

with n = This state is also an eigenstate of the operator Dqq with eigenvalue F{9,K)n 
where 

:r^x Ir^^ siniKT^cosO) . . , 

which has been calculated using ((61). The corresponding transmission spectrum will be 
proportional to the photon number, which is given by 

This has been plotted in Fig|2b]with the angle 9 and detuning Ap/6o. 

Let us first point out that from the left and the right side, the intensity plot is strikingly 
similar to the real space intensity variation in classical light wave diffraction from a straight 



edge 



32 



33| even though intensity variation in these two cases are function of completely 
different set of physical variables. We shall here briefly explain this apparent similarity 
inspite of these differences. Here we have plotted the variation of the photon number as 
a function of the angle 9 and the cavity detuning. Therefore the plot is not an intensity 



plot in the real space. At each value of 6, we obtain a maximum intensity at that value 
of the dispersion shift which corresponds to the number of atoms illuminated in the lattice 
at that angle. As the number of illuminated atoms changes with the change in the relative 
angle 6 so the dispersive shift takes different values depending on the atomic arrangement. 
Nevertheless, the similarity stems from the fact that the factor Dqq which is written in the 
following form, 

1 6 
Doo= 2 + cos{47rjcos^-)]n, (10) 

j=l:K 

mathematically has a similar form of the Fresnel integral, encountered in the intensity profile 
for a straight edge diffraction pattern. There the intensity is a function of C(r) given by 

C(r) = [ cos^Tix^ /2)dx (11) 

where r is dependent on the geometry of the system including the distance from screen. An 
increase in r implies, the evaluation of intensity at a point farther from the straight edge. 
The oscillating behavior of the intensity can be attributed to the functional dependence of 
Eq.l fTTI) . However, it is to be noted, that Eq.l fTOj) involves a summation over the illuminated 
sites K, which is constant and the variation is plotted with respect to 9, which is the angle 
made by the cavity mode with lattice. This summation can be understood in the context of 
the lattice being discrete. But for a fixed K, as we change 9, we are effectively changing the 
number of illuminated sites, as mentioned earlier, and hence the nature of the plot seems 
similar. 

Since the dispersion shift is an indicator of the refractive index of the medium, the 
above result suggests that the refractive index of a given quantum phase is dependent not 
only on the site distribution of the atomic number, but also on the angle between the 
propagation direction and the optical lattice. This is a unique feature of this system. It 
may be recalled that in well known optical phenomenon like Raman Nath scattering due 
to diffraction through a medium with periodically modulated refractive index 3J] or in 
Brillouin scattering in non-linear medium 35| , there is also a frequency shift due to dispersion 
of the transmitted electromagnetic wave through the medium. However the mechanism of 
the dispersion shift as a function of angle between the cavity mode and optical lattice as 
explained in the preceding discussion is fundamentally different from these cases. 
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FIG. 2: (Color Online) (a)Tlie top part of this schematic shows the way atoms are present at the 
intensity maxima regions of the illuminating cavity mode when the optical lattice is illuminated at 
0°. However, changing the angle by 60° results in the decrease of the number of illuminated sites 
to half, (b) Variation of the photon number Q (color axis) with detuning ^ and 9 {in degrees) 
for N=M=30, K=15, k =0.55q, when the atoms are in a MI state and are illuminated by a single 
standing wave cavity mode. 
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2. Superfluid 



Next we consider the case when the atoms are in SF phase. The SF wave function in the 
Fock space basis can be written as superposition state, namely 



1 / jv' 

(rij) 

where nj denotes the number of atoms at the jth site while {nj) denotes a set of nj for a 
particular Fock state. Unlike the MI case, here Dqq acts on a superposition of Fock states 
each of which is an eigenstate of this operator. Each such Fock state carries a different set 
of |ni, 77-2, .., ?t,m)- Hence, 

(^lajaol^) = r o (13) 



mln2l..nMl (Ap - 6oF,{e,K,n,)y + 



where, Fs{6, K,nj) is the eigenvalue of the operator Dqq acting on a particular Fock state. 
Here the Fg functions are generalization of the F function described in ([8]), for the case of 
SF phase in which the number of particles in each site is different as, 

Fs{6,K,nj)= cos'^{jTTCosd)nj (14) 

j=l:K 

where j is the site index and rij is the occupancy of site j. It can be easily checked for 
rij = n for all j, Fg = nF. 

The decomposition of the many body states in Fock space basis is now mapped in the fre- 
quency shifts of the cavity mode. The probabilistic weight factor is mapped in the intensity 
of the peak at that particular value of dispersion shift. At a particular angle of incidence, the 
singular peak of MI now breaks into multiple peaks with varied peak strengths and disper- 
sion shifts. Each particular peak corresponds to a particular group of Fock states which have 
the same value of Fs{6, K,nj) as defined in Eq. (fT4|) . Now as the angle 6 is being changed, 
the effective number of illuminated sites change. This changes the value of Fs{6, K,nj) as 
well as the set of Fock states which yield the same value of Fs{6, K, rij). 

This can be understood clearly by taking a case where the number of Fock states involved 
is small. The corresponding states is few body correlated state which is a few body analogue 
of a superfluid state, where the number of Fock states involved is thermodynamically large. 
Consider the case when 2 atoms are placed in 3 sites among which the first 2 sites are 
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illuminated through a cavity mode. Here the average occupancy per site nj is 2/3. At 
9 = 0°, Fs{9, K,nj) = Ylf=i^j- Hence states yielding the same amount of shift would be 
|2,0,0) , |0,2,0) and |1,1,0) and the corresponding value of Fs{6, K,nj) is 2nj. On the 
other hand the states in which only a single atom is present in the illuminated sites, such as 
|1, 0, 1) or |0, 1, 1) will show a shift of just uj in the cavity frequency. However, |0, 0, 2) does 
not show any shift. Hence, the Fock states get distributed into groups having 3, 2, 1 states 
respectively, each group having a different value of the dispersion shift. 

Now as the angle between the cavity mode and lattice is varied, the effective number of 
illuminated sites change and so changes the set of Fock states that has same Fs{6, K,nj). 
For example, at an angle 60°, only the alternate sites are effectively illuminated. Therefore 
in this case it is just the central site which is illuminated. Now under these circumstances, 
state |0, 2, 0) show a distinctively separate shift from the states |2, 0, 0) and |1, 1, 0). It may 
be recalled that at = 0° all these states had the same shift. Moreover, now the states 
|0, 1, 1) and |1, 1,0) will have same shift since only one atom gets illuminated in this case. 
Thus this evolution of Fock states is imprinted in the shift of the cavity frequency. This is 
shown in FigJ3al The transmitted intensity from the cavity mode will be proportional to 
the photon number in that cavity mode. The color plot depicts this photon number whereas 
the black lines (explained in the legend of Figj3K) correspond to the Fs{6, K,nj)6o for each 
Fock state as a function of 6. As can be inferred wherever there is a maximum overlap 
of Fs{6, K,nj), the same location corresponds to the intensity peak. However an increase 
in the cavity decay rate k shows that the individual peaks cannot be resolved as shown in 
FiglHH 



B. Double mode 



We shall now consider the case where two cavity modes are excited. The corresponding 
photonic annihilation operators are given by clq and ai. Following ^ we also assume that 
the probe is injected only into do, and hence, rji = 0. Also both have the same frequencies, 
i.e., Uq = Ui and are oriented at angles 9o and 9i with respect to the optical lattice. From 
Eq.©, ai = yields 



ai 



irjoSiDioe- 



([Ap - + Vlrn)] + ifi)i[\ - i^J-m - ^m)] + if^) 



(15) 



13 





to 



(a) 



90 

8 

(b) 



180 



FIG. 3: (Color Online) Variation of the photon number (113p (color axis) with detuning -j^ and 
9{m degrees) for N=K=2, M=3, k=0.5(5o when the atoms are in SF state and are illuminated by 
a single standing wave cavity mode. In (a) with individual Fock states, represented by the black 
lines for k =0.15q. In (b) for k =Sq, the spectrum is blurred. 



where 



|(Dn + ^oo) 



00 J 
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are now operators acting on the Fock space. This leads to 



(16) 



a\ai 



(17) 

Here, Dqi = dIq and the expectation value of a[ai gives the photon number in this mode. 

The above problem is equivalent to two linearized coupled harmonic oscillators which 
show mode splittingls^. Briefly two such harmonic oscillators with natural frequencies Ui 
and U2 coupled to each other by a perturbation (, is described by the following set of coupled 
equations. 

dxi . 

— = -lUiXi + CX2 

at 



dx2 



-iU2X2 + (Xi 
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The normal modes of such a system are 



.= ^±^(^P + C (18) 

In the current problem, the shifted frequencies are given by the eigenvalues of 

Urn ± (19) 

acting on a particular state of the system. A particular case of interest will be when 9q = 6i, 
this implies Dqq = Dn = Diq = D. Then the photon number djai is 

a ai = ^ 20) 

[(Ap-2fi^)2 + «:2][A2 + «:2] 

Thus one of the normal mode is independent of the atomic dispersion, however the other 
mode disperses by twice the value for a single mode. We shall now consider the case when 
the cavity modes are Standing Waves. The many body ground state shall be considered as 
either a MI or SF phase. 



1. Mott Insulator 



Again we shall first calculate the two mode transmission spectrum when the cavity con- 
tains atomic ensemble in a MI state given in ((71). The operators Dqq and Du are given by 
([6]) with eigenvalues F{dQ, K)n and -F(6'i, K)n. Also for such a MI state the operator Z)io is 
given by 

cos{jTccos9o)cos{j7icos9i)nj (21) 

j=l:K 

When this operator Diq acts on ([7]) its eigenvalue is given by F{9o, 9i, K)n, where 



F{9o,9i,K) 



sin{KTT 



cosdo+cosdi ' 



sinin 



cosBq+cosB 



—^C0S{{K + 1)7T 



C0s9q + 0056*1 , 



+ 



sin{KTT 



cosOq—cosOi ' 



sm vr 



cosOp—cosd 



—^C0S{{K + 1)7T 



cos9q — cos9i 



The photon number is hence given by 



{^Milalail'^Mi) 



(22) 
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where, 

/ = i^MiloOml'^Mi) = ^-^{F{eo,K) + F{e^,K))n (23) 

^ = i^Mjl^ml^Mi) = n^ ^'^^^'" ; ^^' ^ + SmOo, ^1. ^ (24) 

are the eigenvalues of the operators u and Q acting on MI state respectively. The normal 
modes are hence given by f^J-", and therefore the amount of mode splitting is given by 2 J-". 




(a) (b) 

FIG. 4: (Color Online)(a) shows the variation of \F{6q,6i, K)\'^ (color axis) with the angles 6q and 
6i. (b) shows the mode splitting 27^ Eq. (j24p variation (color axis) in units of di with angles 9q and 
01, N=5, K=5, M= 5. In both cases, 6*0 and Oi varies from 0° to 180°. 

The photon number and the mode splitting, are given by the expressions (12 2 p and fr24p 
respectively. Both are dependent on the value \F{9q,9i, K)]"^ and are thus related to each 
other. Fig. 1 (a) shows the plot of the |F(^o,^i,^)P for n = 1 MI state and Fig. H (b) 
shows the mode splitting at specific values of 9q and 9i and thus very clearly demonstrates 
their inter-dependence. 

This relation is also refiected in the plots of resulting transmission at certain demonstra- 
tive values of 6*0, 9i as plotted in Fig. Oand as explained below. 

First we consider the case when both 9o and 9i are being varied from 0° to 180° always 
maintaining the relation 9o = 9i. The corresponding F[9o,9i, K) function shows a number 
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of maxima along the line = in Fig. HJ^a). The photon number here is given by the 
expressions (120} [22]) where it was seen the normal modes will be zero and . In Fig. |5]|^a) 
we have plotted this variation in photon number along the color function of Q\ and 

Ap/(5o. Therefore at each Q\ one gets a maxima at a value ^ = and at twice the value for 



a single mode case (Fig. |2b 
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(d) 



FIG. 5: (Color Online) Variation of the photon number (I22p (color axis) with detuning and 
the angles and 0i(in degrees) for N=M=30, K=15, k;=0.55o when the atoms are in MI state, 
for double standing mode case. In (a)0o=^i) k=0.5(5o. (b) shows the variation when 0o=O°, Q\ is 
varied and k = b. Similar cases are shown in (c) with ^o=60° and in (d) when Qq= 90° 
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However, in Fig. [5](b), (c) and (d), we describe the case when 6q is kept constant, while 
the other angle 6i is constantly being varied. Corresponding plots show that the maximum 
number of photons scattered from one mode at an angle Oq will be collected by ai only when 
6i = ±6q or tt ± ^^0- When 6i = 6q, the second mode is parallel to the first mode. When 6i 
= -6q, angle of scattering is equal to the angle of reflection. This is also observed at tt ± 
[l^ . It is at this co-ordinate, the \F{6o,6i, K)\'^ , mode splitting as well as the transmitted 
intensity will show maximum behavior. This is also seen from the 6q = 9i and 9q= -k — 9i 
lines in Fig. |H For example in (c),when 9q = 60°, the plot shows maximum mode splitting 
and intensity at 9i =60° and 120°. 

One can also study the diffraction pattern of such system in the limit where the shift in 
the cavity frequency due to dispersion given in f[T7|) is neglected. In that case the transmitted 
intensity will be directly proportional to the eigenvalue of dIqDiq. This particular limit has 



been explored in 



18l . |20| for the two mode case and shown to consist of two parts. The 
first part is due to classical diffraction and second part shows fluctuations from such classical 
pattern. The above analysis in this work suggest an enrichment of these diffraction features 
to a considerable extent once the frequency shift due to diffraction is taken into account. 



2. Superfluid 

Now we consider that the cold atomic condensate is in SF ground state. The transmission 
through the SF can be obtained by taking the expectation value of the photon number 
operator f[T7j) in a SF state. This gives 

1,.^ ^ V 6l[F,{9o,9„K,n,)]'\r^o\' 

(25) 



where 



/„, = |(F.(«i,/s:,%) + f.(»o,A-,n.,)) 



are respectively the eigenvalues of (2jm and Qm operators acting on a particular Fock 
state. These are in terms of Fs{9q, K,nj) functions which were first described in ( IT^ . 
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Fs{9o, 9i, K, Hj) is given by 



Fs{0o,9i, K,nj) = cos{jTTcos9o)cos{j7rcos9i 

j=l:K 



(26) 
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FIG. 6: (Color Online) Variation of the photon number (I25p (color axis) with detuning ^ and 
angles and ^i(in degrees) for N=K=2, M=3, k=0.55o when the atoms are in SF state for double 
standing mode case, with individual Fock states depicted by black lines; in (a)0o=^i (b) Oq= 0° 
{c)9o= 30° {d)9o= 60° 

In Fig. [6] (a), the cavity modes are oriented at the same angle with the lattice axis ie. 
6*0 = 6*1, and are together varied from 0° to 180°. Here Fs{9q, K,nj) = Fs{9i, K,nj) = 
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FsiOo,Oi, K,nj) and thus corresponds to the case described in Eq. (l20|) . In this figure, the 
photon number has been plotted against the angle 6i and Ap/6o- The plot exhibits that at 
each angle there is maxima in the photon number when the value of the dispersive shift is 
either zero or twice its corresponding value for the single mode case. 

FigJH] (b), (c) and (d) describes the variation of the photon number with the angle 6i 
and Ap/6o, when one of the modes, 6q is kept at a constant angle while the other, 6i is 
changed from 0° to 180°. This makes Fs{9q, K, nj), Fs{9i, K, rij) and Fs{6q, 9i, K, rij) change 
separately for individual Fock states. To understand this, we again consider the case when 
2 atoms are placed in 3 sites among which 2 are illuminated. We will have 6 Fock states and 
the average occupancy will be 2/3. However each Fock state now gives maximum intensity 
for two values of Ap/ 6o corresponding to the two normal modes fn. ± J^nj ■ 

FigMjo), shows the intensity distribution when = 0°. Here, we observe that when 
9i=0°, the 6 Fock states distribute themselves into groups having 3,2,1 states and each 
group has a distinct value for the dispersion shift depending on the occupancy. However in 
FigJ^c) and ini^d), where the angle Oq = 30° and 60° respectively, and the angle 9i is varied, 
we observe the Fock states separate out and each takes a distinct value for the dispersive 
shift. This allows us to identify each Fock state separately. In Figj6]^d), it is to be noticed 
that at specific values of 6i, the intensity drops to zero, although there is a superposition of 
the Fock states. This happens as the intensity depending on Fs{9q,9i, K,nj) becomes zero 
at these particular values as can be seen from eqs. ( 12511261) . 



IV. TRAVELLING WAVES 

Fig. [7] depicts our model system where an optical lattice is shown to be illuminated by 
two ring cavities. We have considered that these cavities allow the waves to propagate only 



in one direction. Such cavities generate travelling wave modes 136-381. These modes are 

n n 

described by [8|, |36|, u{rj)Tw = exp(z(fc • Vj + 0)) where is constant phase factor which 
has been set to zero. For such waves, the operator Dqq becomes just 

^00 = u^uirij = ^ fij (27) 

j=l:K j=l:K 

as u'^ui = 1. Thus the eigenvalue of Dqq for a given Fock state will be nK which is just the 
number of atoms in the illuminated sites. Thus the dispersive shift in single mode case will 
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not depend on the angle 9. 




FIG. 7: (Color Online) Schematic diagram of the atom cavity system for travelling wave. The 
optical lattice is created from two counter propagating laser beams and has a site spacing d. The 
two ring wave cavity modes, MODE 1 and MODE 2 are at angles 6q and Oi respectively with the 
axis of the optical lattice. The MODE 1 is being pumped by a pump laser with amplitude ??o while 
MODE 2 is not being pumped but is used to collect the scattered photons by a detector. In the 
single mode case, the detector is also fixed in MODE 1, and has not been shown in this figure. The 
ring cavities, are set in a way that the waves is allowed to propagate only in one direction. 

However this is not the case when two cavity modes are excited. As we have seen for 
the case of standing wave, the mode splitting which in turn influences the transmission 
through such cavity is closely related to the relative angle between the two modes through 
the function F{6o,6i, K) . In the current case the mode splitting is also dependent on the 
relative angle between the cavity modes since only the eigenvalue of the operator Diq is 
angle dependent which is given by. 




(28) 



j=l:K 



21 



A. Mott Insulator 



Again, we first consider the cold atomic condensate in a MI state. The eigenvalue of Cj 
and Ct f[T^ when acting on the MI state ([7j) are 

g = = = nK6i 



g = {^^,\Q\^) = J ^^^^ ''^^' y+\G{9o,9,,K)n5,\^ = \G{9,,9r, K)\n5, (29) 



Here G{9o,9i,K) is 



G{9,M.K)= _,\,,,,„.,,„„L,^/ (30) 



This system is equivalent to two coupled linearized harmonic oscillators fllSp . but with 
same natural frequencies ie., uj\ = uj^ = ujo and coupled by a perturbation C,. The normal 
modes for such a system is given by Wo ± C. In the current problem, the normal modes are 
hence given by (7 ± ^ and therefore the amount of mode splitting is 2Q. 

The photon number (fTTl) is, 

{^Mi\a%\^Mi) = ([A,-(, + ^)P + i?)aA,-(,-^)P + .^) 
Fig. [8] (a) depicts the variation of function G{9q,9i, K = 5) with 9i and 9o. For a 
particular value of 9o and 9i, this function takes the maxima value when the argument of 
the function, ie., {cos9q — cos9i) will become zero, ie., when 9i = ±^o-This can be seen from 
the 9o = 9i line. 

Fig. M (b)-(d) depicts the variation of intensity with Ap/60 and 9i for a fixed value 
of 6*0. The plots show two symmetrically placed transmission peaks, whose separation is 
again proportional to G{9o,9i, K) and therefore will also show a maxima when 9o = ±9i. 
Physically 9i = 9q corresponds to the case when both the ring cavities are oriented at 
the same angle, while ^0 = corresponds to the case when scattering is at the angle 
of reflection. However, in Fig. |8] (b), when 6'o=0°, we observe an additional maxima at 
9i = 180° because, at this value, function G{9o,9i, K) = ^""^f^^-* also shows a maximum 
behavior (Fig. E^a)). Also it is clearly seen that in all these plots, both the normal modes 
symmetrically vary around the average value i.e., nK. This average value is shown by a 
dotted black line in the Fig. [8] (b). As clearly seen, it is independent of the angles between 
the lattice axis and the cavity modes, and is only dependent on the total number of atoms 
present in the illuminated sites. 
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FIG. 8: (Color Online) (a) Variation in G{9(),9i, K = 5)(color axis) with and 0i(in degrees), 
(b) Variation of the photon number (color axis) with detuning ^ and 0i(in degrees), N=M=30, 
K=15, K =0.5(5o when the atoms are in MI state, single travelling mode case. ^0=0° {c)9o= 45° 
(d)0o= 90° 



B. Superfluid 

In this case also, the mode splitting only depends on the eigenvalue of the operator Diq, 
but the eigenvalues are different for different Fock states. Fig. |9] depicts the same for 
travelling mode case. 
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FIG. 9: (Color Online) Variation of the photon number(Color Axis) with detuning and 9i{m 
degrees), N=K=2,M=3, n =0.16o when the atoms are in SF state, two mode case for travelling 
wave cavity with individual Fock states depicted by black lines in {a)6o= 0° (b)6'o= 45° {c)6o= 60° 
{d)9o= 90° 



The photon number for this case will be given by, 



(*|ajo,|*) 



-y 



{Gs?\rio\ 



n^\n2\...nM\ {[\ - {g + QsW + - {d - Qs)f 



(32) 



where g = nKSi and Qs = G{9o,9i, K,nj)6i where, G{6q,9i, K,nj) is the eigenvalue of Diq 
operator on a Fock state with rij particles on the j-th site. It may be again noted that in 
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a SF state nj varies with the site index j for a given Fock state. The transmission spectra 
is shown in the Fig. ([9]) for certain demonstrative values of 60, Gi- In each case we have also 
plotted the Qg for different Fock states which are superposed to form the superfluid. 

V. CONCLUSION 

In our work, we have analyzed cold atomic condensates formed by bosonic atoms in an 
optical lattice at ultra cold temperatures. It has been suggested that such system when 
illuminated by cavity modes, can imprint their characteristics on the transmitted intensity. 
We have studied the off resonant scattering from such correlated systems by varying the 
angles that the cavity modes make with the optical lattice and thus obtained the transmission 
spectrum as a function of the detunings and the dispersive shifts. 

The main result of our work reveals the pattern in the shifts of the cavity mode fre- 
quency as the relative angle between the cavity mode and the optical lattice is changed. 
As we have pointed out in section IIII A II that a change in the dispersion shift implies the 
effective change of the refractive index. Thus our finding implies even for a given quantum 
phase, as the relative angle between the mode propagation vector and the optical lattice 
changes, the cavity induced dispersion shift or the effective refractive index of the medium 
also changes. This highlights the uniqueness of such quantum phase of matter as medium 
of optical dispersion. 

For the single mode case discussed in section IIII At in MI phase, we have seen that the 
transmitted intensity depends on the number of atoms in the illuminated sites, since the 
presence of an atom shifts the cavity resonance and this shift is directly proportional to the 
number of illuminated atoms. The SF phase is however a superposition of many Fock states 
and set of Fock states group correspond to same shift. However changing the angle, these 
group of Fock states change thus providing more information about the system. 

As discussed in next section IIII Bl when two cavity modes are considered, the system 
shows mode splitting between the cavity modes coupled by the atomic ensemble. This was 
clearly visible in the MI case. In the SF state, at some specific angles of illumination, the 
Fock states of SF distinctly map to different frequency shift. Thus giving the Fock state 
structure of the system. However, it was noticed that such a system can only be achieved 
through high finesse cavities, as such characteristic features in the plots for the SF phase 
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become blurred for an increase in ^ values. Such system when illuminated by ring cavities 
show different features of intensity transmission as shown in the section IIVI that describes 
the situation where the cavity modes are travelling waves. Thus the nature of diffraction 
pattern is also dependent on the nature of the quantization of the electromagnetic wave 
inside the cavity, which was also mentioned in the earlier work on such quantum diffraction 



by fermionic cold atoms 31 1. 



As our analysis shows that the variation of the relative angle between the cavity mode 
and the optical lattice can resolve the Fock space structure of a quantum phase in detail, 
by varying the effective number of illuminated sites. It has been pointed out in experiment 



described in ref. 10| that it is possible to study the correlated many body states of few ultra 



cold atoms in such cavity within the currently available technology. A few body correlated 



system of ultra cold fermions was also experimentally achieved recently [39|. For example 
in Fig. [3] and Fig. |6] it has been shown that in the limit of small cavity decay rate k and 
for few number of particles in the illuminated site, in a superfluid phase or more correctly a 
few body analogue of a superfluid state it is possible to identify the extent of superposition 
in different parameter regime. Such identification is potentially helpful in various type of 
many body quantum state preparation. 
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